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Abstract 

Response of pure Ising systems to time-dependent external magnetic 
fields, like pulsed and oscillating fields, are discussed and compared here. 
Because of the two time scales involved, namely the thermodynamic re- 
laxation time of the system and the pulse width or the time period of 
the external field, dynamically broken symmetric phases appear sponta- 
neously when both become comparable. A particularly simple case is that 
of an Ising ferromagnet below its static critical temperature, when it is 
perturbed for a short duration by a pulsed magnetic field competing with 
the existing order in the system. If the field strength and duration is more 
than the threshold (dependent on the temperature), the system, and con- 
sequently the magnetization, switches from one minimum to the other of 
the static free energy. This magnetization reversal transition here shows 
intriguing dynamic transition behaviour, similar to those for oscillating 
fields. Monte Carlo studies for such dynamic transitions are discussed 
and compared with the mean field results for the same and the Monte 
Carlo results for the oscillating field case. In particular, we discuss about 
the Monte Carlo results for the fluctuations and their growth behaviour 
near this magnetization reversal (dynamic) transition point. 

1 Introduction 

The response of pure ferromagnetic Ising systems under the action of time- 
dependent external magnetic fields are being studied intensively these days. A 
whole new class of dynamic phase transitions emerged from the study of such 
driven spin systems under different time dependences of the driving field. Orig- 
inally it was seen that when a pure Ising system is driven by an oscillating 
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external field, it settles to a dynamic steady state with spontaneously broken 
symmetry when the time period of the external field becomes smaller compared 
to the system's relaxation time. In this steady state, the period averaged mag- 
netization acquires a net non-zero value, while the period averaged magnetic 
field remains zero jl], ||, ||, Q| . The universality class of this dynamic transition 
has recently been established, from extensive Monte Carlo studies, to be the 
same as that of the equilibrium transition in the pure Ising system |Q, |). Later, 
efforts were made to investigate and analyse the response of such systems under 
a pulsed magnetic field which, for a short duration, either favours the exist- 
ing magnetic order of the system (positive pulse) or competes with it (negative 
pulse) . The results of the positive pulse case were analyzed by extending appro- 
priately the finite size scaling technique to this finite time-window case, and it 
did not indicate any new phase transition or thermodynamic scale ||. However, 
the negative pulse case indicated a new dynamic transition and thermodynamic 
scales, as the appropriate pulse width and strength just manages to help the 
system overcome the "double-well free-energy barrier" and the magnetization 
reversal transition occurs |?|, [| [1. 

In the next section we briefly introduce the model and the magnetization 
reversal transition in this negative pulse case. We also briefly summarize here 
a few analytic results obtained earlier in the mean field limit. In section 3, we 
report Monte Carlo results for the fluctuations near this magnetization reversal 
transition in two dimensional Ising systems. In the concluding section, we anal- 
yse these observations and also compare them with those obtained for the Ising 
model under oscillating field. 

2 Model 

The model we consider here is the Ising model with nearest-neighbour interac- 
tion under a time-dependent external magnetic field, described by the Hamilto- 
nian 

M « 

Here, J denotes the strength of the cooperative interaction between nearest 
neighbour spins put under a time-dependent field h{t). Si and Sj denote the 
spins at lattice sites i and j respectively and the nearest-neighbour spin pairs 
are denoted by [...]. We consider such a system at temperatures T below its 
static critical temperature T c . The external field is applied at t = to after 
the system is brought to equilibrium, having an equilibrium magnetization mo 
(T). The external field h(t) has a time dependence given by h{t) = —h p for 
to < t < to + At and h{t) = outside this range. Typical responses of the 
time-dependent magnetization m(t) of the system under different h(t) are in- 
dicated in Fig. 1. For appropriate combinations of h p and At, at a fixed T, a 
magnetization reversal transition occurs when the magnetization of the system 
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eventually switches to —mo. It may noted that a magnetization reversal phe- 
nomenon trivially occurs in the limit At — ► oo for any non- vanishing value of h p 
at any temperature T < T c . This is however a limiting case of the transition, 
which is studied here for any finite value of At . Also, in the case of magneti- 
zation reversal transition, the magnetization reversal need not occur during the 
presence of the external pulsed field. In fact, the time required for the system 
to reach its final equilibrium state diverges as the transition point is approached 

0. 

A qualitative understanding for these diverging time and also length scales 
at this magnetization reversal transition point can be easily understood from a 
mean field study of this problem |Q, Using mean field approximation for the 
dynamics, we can write the equation of motion for the average magnetization 

rrii(=< Si >) as 

^i = - mi + tmh ( J S."f "» ), (2) 

where < ... > denotes the thermodynamic (ensemble) average. This equation 
may be linearized near the magnetization reversal transition point and for uni- 
form magnetization one gets 



m(t) =m - - mo 



f AT . 
exp < — [t - t ) 



(3) 



as a solution of Eq. (§). Here AT = T c mf - T, where T c mf = J(q = 0) is the 
critical temperature in the mean field approximation and J(q) is the Fourier 
transform of the interaction J. From the solution (||) the value of mo is always 
restricted to be less than h p / AT in order to make the linear approximation valid. 
Hence during the presence of the field, m(t) will decrease in its magnitude from 
m(io) = too to m(to + At) = m w at the time of withdrawal of the pulse. Since 
there is no fluctuation here, the magnetization will relax back to its original 
value too if m w is positive or relax to a value —too if m w is negative. This can 
be easily seen by solving eq. (||) for t > to + At where h(t) — 0, using m w as 
the initial value of the magnetization. In this time region one gets a Debye type 
relaxation of magnetization to its final equilibrium value, with a relaxation time 
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It diverges at the magnetization reversal transition point, where m w vanishes. 
The prefactor, giving the divergence of r at the static mean field transition 
temperature, is responsible for the usual critical slowing down phenomenon at 
the static transition point (h — 0). The other factor gives the diverging time 
scale of interest here, at any temperature below the static transition tempera- 
ture, where magnetization reversal occurs or m w vanishes due to the appropriate 
combination of the external field pulse strength and duration. A similar solution 
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for the Fourier transform of the linearized version of the equation (g) , with an 
additional derivative with respect to the field h for the susceptibility x(q)> gives 



Here also, the prefactor in £ gives the usual divergence at T™ f , while the other 
factor gives the divergence of interest here at the magnetization reversal tran- 
sition point. Such divergences of r and £ across the dynamic transition phase 
boundary which gives the optimal combination of h p and At at any particular 
T, are also seen || in the Monte Carlo simulations. 

We believe, a similar transition and the associated length and time scales are 
responsible in an essential way for the spontaneous dynamic symmetry breaking 
transition in the Ising systems under oscillating fields |Q, ||. The growing time 
scale induced by the periodically changing field competing with the existing 
order in the system at that point, locks the system dynamically in a broken 
symmetric phase. In order to compare these scenario, we have studied here 
the growing fluctuations in the order parameter m w and the internal energy 
E, as one approaches the magnetization reversal phase boundary in a Monte 
Carlo study. We find that the fluctuations in both these quantities diverge very 
prominently at the magnetization reversal phase boundary. We also indicate the 
behaviour of the derivatives (dm w /dh p ) and (dE/dT), which show some peak 
across the same phase boundary. We compare our results with those reported 
for the dynamic transition under oscillating field 0). 

3 Monte Carlo study 

We study here the behaviour of the fluctuations near the magnetization reversal 
transition in a ferromagnetic Ising system on square lattice of size L, under 
negative or competing field pulses, using Monte Carlo simulations with single 
spin flip Glauber dynamics. Working at temperatures T below the static critical 
temperature T c (~ 2.27 in units of J ), the system is brought to the equilibrium 
state by evolving the system without any external field, from an initial config- 
uration (say all spins up). The equilibrium relaxation time required for this 
system to evolve to its equilibrium state with magnetization uiq is noted. The 
field pulse of strength —h p for a duration At (measured in Monte Carlo steps 
or MCS), competing with this order, is applied at t = to only after the system 
has reached its equilibrium. The magnetization m(t) starts decreasing (from 
its equilibrium value m ) due to the effect of the competing field — h p during 
the time period At. The average value of the magnetization m w at the time 
of the withdrawal of the pulse is noted. The phase boundary of this dynamic 
transition gives the appropriate combination of h p and At that induces the mag- 
netization reversal by bringing m{t) from m(io) = mo to m(to + At) = m w = 0. 
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Therefore across the phase boundary m w changes its sign JtJ], Obviously the 
phase boundary changes with T. We study the behaviour of various thermody- 
namic quantities across these phase boundaries. The thermodynamic quantities 
are averaged over different initial configurations of the system. We note the 
fluctuations over these average quantities as well. In our simulations we have 
considered three different system sizes, viz., L =256, 512 and 1024 for all our 
studies reported here. The typical number of averages is 200 for L = 1024 and 
1000 for L = 256. 

In Fig. 2 we show the growth of fluctuations Am w = L 2 (< > — < 
m w > 2 ) in m w across the phase boundary at two different points, both driven 
by the pulse strength h p . The insets here show the corresponding variations of 
the pulse- withdrawal-time magnetization m w , which crosses zero at the phase 
boundary. Figure 3 shows the growth of fluctuations AE to t — L 2 {< E 2 ot > — < 
Etot > 2 ) in the average total energy E tot per spin pair (in units of J) across 
the phase boundary, again at two points, both driven by T. The insets here 
again show the corresponding m w variations. Both these fluctuations Am w and 
A£" t ot clearly diverge at the phase boundary, in the large system size limit. Fig. 
4 shows the variation of the fluctuation AE coop = L 2 (< E 2 oop > — < E coop > 2 ) 
in the cooperative part E coop (excluding the field contribution) of the total 
energy E to t for the same parameter values as those of Fig. 3(a). Although the 
fluctuation in E coop also grows in general around the transition point, a sharp 
dip is observed there. This dip is identified to be due to the sharp growth of 
-Scoop at the transition point, as shown in the inset. 

In Fig. 5, we show the behaviour of the derivatives dm w /dh p (in (a)) and 
di?tot/dT (in (b)) across the phase boundary at two different points. Unlike in 
the case of the dynamic transition under oscillating fields []l(J , we do not observe 
any divergence in these quantities at the transition point. This is because the 
curves are independent of system size and hence the peak values do not grow 
with L. However, the peaks become sharper at lower temperatures (higher 
values of h p ) as the system approaches the tricritical point occurring on the 
phase boundary 

4 Concluding remarks 

We have studied here the growth of fluctuations and their system size depen- 
dence across the dynamic magnetization reversal phase boundary. This tran- 
sition can be driven by the field pulse strength h p or its duration At or even 
by varying T in the presence of a fixed field pulse but always keeping it below 
T c . 'Critical slowing down' and the divergence of length scales etc at this dy- 
namic phase boundary have already been demonstrated earlier ||], using both 
mean field approximation and the Monte Carlo simulations in two and three 
dimensions. Here we have shown the Monte Carlo results for the growth of 
fluctuations, both in the field- withdrawal-time magnetization m w and in the 
total energy .Etot > across the dynamic phase transition boundary in two dimen- 
sion. The growth of these fluctuations with the system size L clearly indicates 



5 



a divergence at the phase boundary. It may also be noted that although we 
have shown the fluctuations for one specific driving parameter, e.g. Am w with 
h p and AE to t with T, similar fluctuations in both the quantities are seen for 
other driving parameters as well. All these clearly indicate a genuine dynamic 
phase transition at the magnetization reversal point of an Ising system under a 
competing field pulse. 
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Figure Captions: 



Fig. 1 Typical time variation of the response magnetizations m(t) for two 
different field pulses h(t) with different strengths h p , but same width At and at 
same temperature T, are shown. The quantities of interest to characterize the 
response magnetization are also indicated. 

Fig. 2 Growth of fluctuations Am w near the transition, driven by h p , are 
shown for two different points on the phase boundary, (a) for T = 1.0 and 
At = 12 MCS, and (b) for T = 1.50 and At = 2 MCS. The insets show the 
corresponding m w variations; m w = at the transition point. 

Fig. 3 Growth of fluctuations AE tot near the transition, driven by T , are 
shown for two different points on the phase boundary, (a) for h p = 1.75 and 
At = 10 MCS, and (b) for h p = 1.25 and At = 10 MCS. The insets show the 
corresponding m w variations; m w = at the transition point. 

Fig. 4 Growth of fluctuations Ai? coop near the transition point at h p = 1.75 
and At =10 MCS. The insets show the (rather sharp) variations of < -E C oo P > 
and < i?coo P > & t the transition point, giving rise to the dip there in AE coop . 

Fig. 5 (a) The variation of dm w /dh p across two different points on the phase 
boundary: for T =1.50 and At = 2 MCS (open symbol) and T = 1.00 and At 
= 12 MCS (solid symbol); (b) The variation of d-Etot/dT across two different 
points on the phase boundary: for h p = 1.25 and At =10 MCS (open symbol) 
and h p = 1.75 and At = 10 (solid symbol). 
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